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. ABSTRACT 

' Decoherence is the main obstacle to the reahzation of quantum computers. Until recently it was thought that 
quantum error correcting codes are the only complete solution to the decoherence problem. Here we present 
an alternative that is based on a combination of a decoherence-free subspace encoding and the application of 
i strong and fast pulses: "encoded recoupling and decoupling" (ERD). This alternative has the advantage of lower 
^ ' encoding overhead (as few as two physical qubits per logical qubit sufRce), and direct application to a number 
of promising proposals for the experimental realization of quantum computers. 

(N : 

1. INTRODUCTION 

^ ' In the quest to construct large-scale quantum information processors, in particular quantum computers, deco- 
QQ , herence is still the main obstacle to realization. Decoherence is the degradation of quantum information due to 
0^ ' inevitable interactions with the environment. Early skepticism^' concerning the viability of quantum compu- 
I tation (QC) in the presence of decoherence was overcome by the discovery of quantum error correcting codes 
■ (QECCs),^"^ that flourished into a comprehensive theory that incorporates all elements of quantum computa- 
\ tion.'^'^ A QECC relies on an encoding of quantum information into the state of several quantum bits (qubits), 
' and a closed-loop, active error diagnosis and correction procedure. In principle, it is possible to correct ar- 
I bitrary errors using sufficiently large QECCs.^ In practice, however, this may require a very large overhead, 

^ • especially in terms of qubit resources.^ This is troubling in light of the substantial difficulties associated with 

^ 1^ , generating and controlling systems with very large numbers of qubits. In addition, the theory of QECCs is rather 
' abstract, in that it presumes that one can execute certain logical operations, but does not refer to the underlying 
I Hamiltonians governing specific physical systems. Hence it is of interest to explore alternatives to QECCs, that 

^ are more economical in qubit resources, and that are direcly tailored to specific quantum computer proposals. 

' Here we review our recent progress in developing such an alternative, "encoded recoupling and decoupling" 

\ (ERD),^^"^^ and report on some new results. ERD is based on a combination of encoding quantum information 
^ ' into a decoherence-free subspace (DFS)^^"^^ (for a review see Ref. 26), and the apphcation of fast and strong 
^ I dynamical-decoupling (or "bang-bang", BB) pulses. ^''"^^ The DFS encoding provides a first layer of 
■ - - ' protection against decoherence, while the BB pulses are used to efficiently reduce the remaining decoherence. The 
encoded "recoupling" part of ERD refers to the application of control operations that enact universal quantum 
logic, in a manner that is fully compatible with the DFS encoding and the BB pulses, and takes into account the 
experimentally available control resources, such as the underlying system Hamiltonian.^^ The utility of ERD 
as a general method for quantum simulation, universal QC, and decoherence suppression has also been stressed 
and explored by Viola. ■^^ Our work builds in part on earlier efforts to combine universal QC with DFS encoding 
(without BB pulses).32.35,46,47 

The structure of the paper is as follows. We begin in Section 2 with a brief formal summary of decoherence. We 
then give in Section 3 a simple example of a DFS encoding, protecting against collective dephasing. Section 4 
shows how to perform universal QC on this encoding. We then briefly and formally review the dynamical 
decoupling method, in Section 5. The methods are then combined in Sections 6,7, where we show how BB pulses 
can eliminate decoherence sources beyond collective dephasing. In Section 8 we present a new and somewhat 
surprising result: the elimination of off-resonant transitions, that will typically be induced by BB pulses, via BB 
pulses. Section 9 shows how to combine all the pieces together, by including logic gates with the BB pulses. 
Concluding remarks are presented in section 10. 



1 



2. DECOHERENCE 

The dynamics of an open quantum system coupled to a bath is formally obtained from the time-ordered evolution 

U{t) =Texp{-i J H{t')dt') (1) 

under the combined system-bath Hamiltonian 

H = {Hs + He) (S>Ib + Is®Hb+ Hsb 
HsB = Y.S^(^B^, (2) 

7 

where / is the identity operator, Hs [Hb) is the internal Hamiltonian for the system (bath) alone, He is an 
externally applied control Hamiltonian, Hsb is the system-bath interaction Hamiltonian, and the {Bj) are 
operators acting on the system (bath). The S*^ play the role of error generators in QC. Making the standard 
assumption of initially decoupled system and bath, one traces over the bath degrees of freedom in order to obtain 
the time-evolved reduced system density matrix: 

Pit) = TrB[U{t) (p(0) pb(0)) UHt)], (3) 

where p(0) [pb(0)] is the initial density matrix of the system [bath]. As is well-known, p{t) is in general a mixed 
state: Tr[p(t)^] < 1, indicating that the environment has decohered the system (see, e.g., Ref. 48). 

3. SIMPLE EXAMPLE OF DECOHERENCE-FREE SUBSPACES 

We briefly review the simplest example of a DFS. Suppose that a system of K qubits (two-level systems) is 
coupled to a bath in a symmetric way, and undergoes a dephasing process. This can be described by the 
system-bath Hamiltonian 

K 

H^on-deph^j2a]®B' (4) 

i=i 

where cr| is the Pauli-z matrix acting on the jth qubit and is an arbitrary bath operator. Under this 
interaction qubit j undergoes the transformation 

|0),--|0), |l),-.e^^|l)„ (5) 

which (after Tr^) puts a random phase (p between the basis states |0) and |1) (eigenstates of with respective 
eigenvalues +1 and —1). This can also be described by the matrix Rz{(j)) — diag (l, e*"^) acting on the {|0), |1)} 
basis. As implied by (4), we assume that the phase has no space (j) dependence, i.e., the dephasing process is 
invariant under qubit permutations. Since the errors can be expressed in terms of the single Pauli spin matrix cr^ 
of the two-level system, this example is referred to as "collective phase damping" , or "weak collective dccoher- 
ence".'*^ The more general situation when errors involving all three Pauli matrices are present, i.e., dissipation 
and dephasing, is referred to as "strong collective decoherence" or just "collective decoherence".^^'^'*'^® With- 
out encoding a qubit initially in an arbitrary pure state — a\0)j + b\l)j will decohere. This can be seen by 
calculating its density matrix as an average over all possible values of (j), 

/■QC 
R,{4>)\^)MRl{<t>)p{<t>)d(l>, (6) 
-oo 

where p{(j)) is a probability density, and we assume the initial state of all qubits to be a product state. For a 
Gaussian distribution, p(0) = (47ra)~^^^ exp(— (j!)^/4a), it is simple to check that 

o-( I'^l' ah*e-\ 
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The decay of the off-diagonal elements in the computational basis is a signature of decoherence. 

Let us now consider what happens in the two-qubit Hilbert space. The four basis states undergo the trans- 
formation 

|0)iO|0)2 ^ |0)iO|0)2 (8) 
|0)i«)|l)2 ^ e'^|0)i«)|l)2 (9) 
|l)i®|0)2 ^ e*^|l)i®|0)2 (10) 
|l)i®|l)2 ^ e2'*|l)i® |1)2. (11) 

Observe that the basis states |0)i |1)2 and |l)i ® |0)2 acquire the same phase. This suggests that a simple 
encoding trick can solve the decoherence problem. Let us define encoded states by |0l) = |0)i ig) |1)2 = |01) and 
= |10). Then the state jf/'i) = a\OL) + b\lL) evolves under the dephasing process as 

IVi) ^ a|0)i ® e^'^|l)2 + be"f'\l)i ® |0)2 = c'^I^Pl), (12) 

and the overall phase thus acquired is clearly unimportant. This means that the 2-dimensional subspace 
DFS2{0) — Span{|01), |10)} of the 4-dimensional Hilbert space of two qubits is decoherence- free (DF). The 
subspaces DFS2{2) = Span{|00)} and DF5'2(— 2) = Span{|ll)} are also (trivially) DF, since they each acquire 
a global phase as well, 1 and e^*''^ respectively. Since the phases acquired by the different subspaces differ, there 
is decoherence between the subspaces. 

For K = 3 qubits a similar calculation reveals that the subspaces DFS3{1) = Span{|001), |010), |100)} and 
DFSsi-l) = Span{|011),|101),|110)} are DF, as well the (trivial) subspaces DFS3{3) = Span{|000)} and 
DFSsis) = Span{|lll)}. 

More generally, let 

Xk = number of O's minus the number of I's (13) 

in a computational basis state (i.e., a bitstring) over K qubits. Then it is easy to check that any subspace 

spanned by states with constant Xk is DF, and can be denoted DFSk{Xk) in accordance with the notation 
above. The dimensions of these subspaces are given by the binomial coefRcients: d = dim.[DFSK{XK)\ = {xk) 
and they each encode log2 d qubits. 

The encoding for the "collective phase damping" model discussed here has been tested experimentally. The 
first-ever experimental implementation of DFSs used the DFS2{0) subspace to protect against artifially induced 
decoherence in a linear optics setting.'*^ The same encoding was subsequently used to alleviate the problem of 
external fluctuating magnetic flelds in an ion trap quantum computing experiment, ^'^ and figures prominently 
in theoretical constructions of encoded, universal QC.^'^'^^'^^ 

4. ENCODED UNIVERSAL LOGIC GATES 

Next we show how to perform universal QC on the encoding DFS2{0) = Span{|01), |10)}. 

4.1. Single encoded-qubit gates 

Let Xi, Yi,Zi denote the standard Pauli matrices af,af,a?, acting on the ith physical qubit (we will use both 
notations interchangeably). In Ref. 12 it was shown that for the code {|0i) = |01),|li:,) = |10)} the encoded 
logical operations (involving the first two physical qubits) are 

Xi2 = -^{XiX2 + YiY2), 

Yu = -^{YiX2 — ^1^2), 

Z12 = ^{Z,-Z2). (14) 
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These operations form an su{2) algebra (i.e., we think of them as Hamiltonians rather than unitary operators). 
We use a bar to denote logical operations on the encoded qubits. In Refs.^^"^^ these logical operations were 
denoted by T{^,, a G {x,y,z}, and a detailed analysis was given on how to use typical solid-state Hamiltonians 
(Heisenbcrg, XXZ, and XY models, defined below) to implement quantum logic operations using this DFS 
encoding. E.g., the term X1X2 + Y1Y2 is the spin-spin interaction in the XY model, and Zi — Z2 represents 
a Zccman splitting. A static Zccman splitting and a controllable XY interaction can be used to generate a 
universal set of logic gates. Similar conclusions hold when the XY interaction is replaced by a Heiscnberg-'^^' 
or XXZ interaction,^'* or even by a Hciscnberg interaction that includes an anistropic spin-orbit term.*^ We 
remark that, as first shown in Rcfs. 46, 47, the various types of exchange interactions can be made universal also 
without any single-qubit terms (such as a Zeeman splitting), by encoding into three or more qubits,''^' ^'^"^^ a 
result that has been termed "encoded universality".^^ 

The two-qubit gate can be expressed as follows: 

Uij{e,(t>i,4>j) = exp{ieX^.X^^) 

= cos6IiIj +ism9X^^ X^. , (15) 

where 

X^ = X cos (f> + Y sin (I), (16) 

and / is the identity operator. In the context of trapped-ion QC,^* the phase (f)i is the phase of the driving laser 
at the ith qubit, while 6 is proportional to the Rabi frequency, and can be set over a wide range of values. 
Introducing the operators 

X^, ^ \{X,X, - Y,Y,), Y, ^ \{Y,X, + X^Y^) (17) 
(denoted i?^-, i?^- respectively in Refs. 12-14) we can express 



Uij{0,<pi,(j)j) = cosOI + ism6{cos AcpijXij + sin A(pijYij + coti^ijXij + iiin^ijYij), (18) 

where =4>i+4>j. It is simple to check that andl^j annihilate the code subspace {|0l) = |01),|1l) = |10)} 
and have non-trivial action (as encoded X and Y) on the orthogonal subspace {|00),|11)}. Therefore, upon 
restriction to the DFS we can write: 

DFS — 

Uij{9,(t)i,4>j) ^ Uij{6,A(t)ij) 

= exp{i9X Aci>ij ) = cos 91 + i sin 9XAci,ij ■ 

(19) 

The fact that Uij depends only on the relative phase A(pij is crucial in the trapped-ion context: this quantity 
can be controlled by adjusting the angle between the driving laser and the interatomic axis, as well as by small 
adjustments of the trap voltages (which, in turn, control the trap oscillation frequency, and hence the qubit 
spacing), whereas it is much harder to control the absolute phase and hence also This is why the 

code subspace {|01), |10)} enjoys a preferred status over the subspace {|00), |11)}. 

The Hamiltonians in Eq. (14) generate the logic operations (19): 

exp{i9Xi2) = Ui2{9,cp,(l))^Ui2{9,Q) 
exp(i0yi2) = C/i2(0,0,(/.+ |) = f7i2(^,|) 

exp(i^Zi2) = exp(i^yi2)exp(i^Xi2)exp(-i^Fi2) 

= f7i2(^,7r/2)C/i2(^,0)[/i2(-J,7r/2). (20) 

The third line follows from the elementary operator identity 

X4, = Xcos(l) + Ysin<j) = e-'<>>^/'^Xe''^^/'^ (21) 

which holds for any su{2) angular momentum set {X, Y, Z}, i.e., operators that satisfy the commutation relation 
[X, Y] = HZ (and cyclic permutations thereof), in particular also the encoded operators {X, Y, Z}. 
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4.2. Entangling gate between pairs of encoded qubits 

It is a well-known requirement of QC that in order to enact a universal set of logic gates (that allow any unitary 
transformation to be implemented) it is suffient to be able to implement all single-qubit operations (as in the 
previous subsection), and to entangle pairs of qubits. In this subsection we discuss the latter requirement. 

In Ref. 61 the following unitary gate was introduced, which is particularly suitable for QC using trapped 
ions: 

= {hhhh - iX^^ X^^ X^^ X^^ ) 

71 (^12^34 - iX i^^^^X ^^^^) 
= eM-i^XA^^Y^^J. (22) 

This gate can be used to entangle two DFS-qubits. It involves simultaneous control over two phase differences 
A(/)j2, A(/)34, and thus control over the motion of two pairs of qubits. The case A(j)^^ = A(j)34 = was used in 
Ref. 60 to demonstrate entanglement of four trapped-qubit qubits, but this choice is not unique. 

For general exchange Hamiltonians a different method is required. We review the encoded recoupling method 
introduced in Ref. 13. The exchange interaction quite generally has the form 



^^ex = E (23) 

cic—x.,y,z i<Cj 



where 

1 . . . ,.^ 1 



(25) 



andJ^=J-±J^. 

The isotropic (Heisenberg) case corresponds to = Jij . The XY model is the case J?- = Jfj , J^?- = . The 
XX Z model is the case J[j = ^Jij 7^ Jfj- A summary of QC proposals that fall into each category can be found 
in Ref. 13. The free Hamiltonian is 

Ho = Y^leiat, (26) 

i 

where is the single-particle spectrum. In each instance of 7?ex (Heisenberg, XY, XXZ) one typically has 
control over only one type of parameter out of the set { J^", e^}. 

Let = J2m-i,2m i*^ ~ -2) ^t), and = (£2^-1 i e2m) /2. Let A and B be two angular momentum 
operators satisfying su{2) commutation relations. Then if follows from Eq. (21) that the operation of "conjugating 
by A", 

Ca ° exp{iB) = exp{—iAn/2) exp{iB) exp{iA7r/2) 

= exp(-iB) (27) 

causes B's sign to be flipped. 

Consider now an XX Z-typc Hamiltonian where the J+ parameters are controllable but and are fixed. 
For simplicity let us consider just the case J~ = 0. Then we can rewrite 

N/2 
m=l 
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where we have omitted a constant term. The important point is now that and satisfy su(2) commutation 
relations. Therefore by using recoupHng through "conjugation by T^" we can selectively turn on and off the 
single-encoded-qubit rotation T^^ and the encoded-Ising interaction Tf^^T^^^. This example of "encoded selective 
recoupling" establishes that encoded universal computation in the XX Z model can be done using control over 
the J+ parameters alone. 

Next, consider the XY model, with controllable J^, but fixed Si. To implement encoded single-qubit op- 
erations, we can use the same encoded recoupling method as for the XXZ model. As for encoded two-qubit 
operations, we now no longer have the (jfcrj terms. Since the XY model with nearest-neighbor interactions can 
be shown not to be universal, we turn on also next-nearest neighbor terms (these can still be nearest- 
neighbor in a 2D hexagonal geometry). First note that Ct^^ o Tlj ~ icrj^CTfrfg. Now assume we can control 
Ji^; then, using conjugation by 7r/4: Cirp,:^ o (Ct^^ o = o-|(c7| — erf )/2. Since uf cr| is constant on the code 
subspace it can be ignored. On the other hand, (Tfcrg again acts as — TfTI , i.e., as an encoded (g) a^. This 
establishes universal encoded computation in the XY model. 

Taken together, the results in this section show how universal QC can be implemented using qubits encoded 
into a DFS offering protection against collective dephasing, while using only reasonable models of physically 
controllable Hamiltonians. We now move on to a discussion of how to reduce additional sources of decoherence. 

5. DYNAMICAL DECOUPLING PULSES 

We briefly review the dynamical decoupling method, introduced in Ref. 27, and further developed in Refs. 15, 18, 
19, 21, 28-45. A set of "symmetrization" operations is chosen such that they form a discrete subgroup of the full 
unitary group of operations on the Hilbert space of the system. Denote this group Q ~ {^}jJi where = I. 
The cycle time is Tc = \Q\At, where \G\ is the number of symmetrization operations, and At is the time that the 
system evolves freely between operations under Uo{t) = exp{—iHst). The symmetrized evolution after a single 
cycle is given by 

151 

U{T,) = II vJUo{At)Vj = e'^««^% 

where the evolution under Hsb + Hb has been neglected during application of the pulses Vj, which quantifies the 
sense in which these pulses have to be fast and strong (for details see, c.g, Refs. 21,27,30,44,45). iJeS denotes 
the resulting effective Hamiltonian, that can be computed to any order using the Magnus expansion. In the 
BB limit one is interested in the evolution U{NTc) after N cycles, such that N ^ 00 and At ^ while NAt is 
finite. As shown in Ref. 34 one can then retain only the first order Magnus term, approximating i?eff by 

1 

^eff = E ^I^v^ ^ ne(if). (29) 

I I j=i 

The map Ug is the projector into the centralizer Z of Q, defined as 

z{g) = {x\ vjxvj = X, vvj e g}. 

Since the pulses gj are unitary the centralizer equals the commutant in this case: 

z{g) = {x\ [x,v,] = o, WjGg}. 

It is clear that Ug commutes with the adjoint action Vj ■ Vj for all j, so that if g is generated by {/, Hs, S-y}, 
the evolution will proceed without the operators affecting the system as errors: the evolution has been 
symmetrized with respect to g.^^ For a geometric interpretation of this symmetrization see Ref. 41, where it was 
also pointed out that while the BB pulses should have the effect of rotating H to the vertices of a symmetric 
object, they need in fact not form a group. 
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6. CREATING COLLECTIVE DEPHASING CONDITIONS USING DECOUPLING 
PULSES: REDUCING DECOHERENCE DURING STORAGE 

One of the important advantages of the DFS encoding {|01), |10)} is that it is immune to collective dephasing. 
However, other sources of decoherence inevitably remain. In this and the following section, we algebraically 
classify all additional decoherence effects and show how they can be eliminated. 

6.1. Creating collective dephasing on a pair of qubits 

First, let us analyze the effect of breaking the collective dephasing symmetry, by considering a system-bath 
interaction of the form 

^d^ph(2) =Zi(3Bt + Z2(3BI (30) 

where arc arbitrary bath operators [compare to Eq. (4)]. This describes a general dephasing interaction 

on two qubits. The source of such dephasing during storage can be long wavelength, randomly fluctuating 
ambient magnetic fields, that randomly shift the relative phase between the qubit |0) and |1) states through 
the Zeeman effect. The interaction can be rewritten as a sum over a collective dephasing term Zi + Z2 and 
another, differential dephasing term Zi — Z2, that is responsible for errors on the DFS: 

Here B^^^ = {Bf + -B|) /2 and B^jj = {Bf — -Bf) /2. If S^jj were zero then there would only be collective 
dephasing and the DFS encoding would offer perfect protection. However, in general B^^f 0, and the DFS 
encoding will not suffice to offer complete protection. 

The crucial observation is that, since Z-^ — Z2 oc Z\2 [recall Eq. (14)], the offending term causes logical errors 
on the DFS.^^ As shown in Refs. 15,37, then the problem of B^^f 7^ can be solved using a series of pulses 
that symmetrize H'^'j^^'"^^ such that only the collective term remains. To do so note that since the offending 
term cx Z12, it anticommutes with X12 = ^(^1X2 + Y1Y2). At the same time X12 commutes with Zi + Z2. 
This allows us to flip the sign of the offending term by using a pair of ±7r/2 pulses in X12, while leaving only 
the collective term. Evolution with the flipped sign followed by unaltered evolution leads to cancellation of the 
offending term. Specifically:^^ 



or, in terms of gates: 



e-'-f^^«^C7i2(-|,0)e-'-^^«^C7i2(|,0) = g-i(Zi+22)®BJ,i2r^ ^33) 

where tJij{9, A(pij) was defined in Eq. (19). This equation means that the system-bath coupling effectively looks 
like collective dephasing at the end of the pulse sequence. Thus, the system is periodically (every 2r) projected 
into the DFS. 

In order for the the procedure described in Eq. (33) to work, the gate C7i2(±|,0) must be executed at a 
timescale faster than the cutoff frequency associated with the fluctuating (magnetic) flelds causing the differential 
dephasing term in H^^'^^''^^ . 

6.2. Creating collective dephasing on a block of four qubits 

So far we have discussed creation of collective dephasing conditions on a single DFS qubit. However, it is essential 
for the reliable execution of an entangling logic gate to have collective dephasing over all four qubits participating 
in the gate, even if only two are coupled at a time. A procedure for creating collective decoherence conditions 
over blocks of 3, 4, 6 and 8 qubits was given in Ref. 15. Here we show how to do the same for a block of 4 qubits 
with collective dephasing. 



7 



Let us start with a general dephasing Hamiltonian on N qubits, and rewrite it in terms of nearest-neighbor 
sums and differences: 



N 

i=l 

N/2 

= (^2j + Z2j-l) ® B+ + {Z2j - Z2j-l) ® B^j, 



(34) 
(35) 



where -B^- = (-62^" ± 52j-i)/2. As noted above, Z2j — -^2^-1 -^2^-1, 2j, so that to chminatc all nearest-neighbor 



differences of the form (Z2j — ■^2^-1) we can use the collective decoupling pulse X„„ = e* 2 ^^^-^-^^ : 



(36) 



or, in gate terms: 



-iHsBT 



N/2 



C/2i-i,2i(-2,0) 



N/2 



^U2j-l,2j{^,0) 



The next step is to eliminate next-nearest neighbor differential terms. To this end let us rewrite the outcome of 
the X^n pulse in terms of sums and differences over blocks of four ions: 

N/2 N/2 

Y,{Z2j+Z2j-l)®B+ = Y.^Z2,+2 + Z2, + l+Z2j+Z2j-^]®B+: + 



N/2 

+ X] [(^2j+2 - Z2j) + {Z2j+1 - Z2J-1) 



B. 



2i 



(37) 



where -825'^ = {B2j+2 ^ ^tj)/^- "^^^ term in the first line contains only the desired block-collective dephasing 
over 4 ions. The term in the second line contains undcsircd differential dephasing terms that we wish to eliminate. 
But these terms once again have the appearance of encoded Z operators, between next-nearest neighbor ion pairs. 
Therefore we need to apply a second collective pulse X„„„ = 6*^^2^-1,2^+1 gif X2j,23+2 ^ ^j^^^ applies encoded 

X operators on these qubit pairs. At this point we are left just with collective dephasing terms on blocks of 4 
qubits, as required: 



,-i2r-Z,^liZ2j+Z2. 



E "ii' {Z2, +Z2,-i)®Bi 



N/2 




)U2J 






'n/2 






i 


(2) 







L(-|,0)t/2,-2,+2(-|,0) 



- TT - TT 

i^2j-l,2j+l( 2 ' 0)f^2j,2j+2( 2 ' 0) 



= e,-^'^rY.'^ll(Z2,+2 + Z2, + i+Z2,+Z2,-i)!!i,B+:+ _ 



(38) 

This pulse sequence is important to ensure that collective dephasing conditions will prevail during the execution 
of logic gates between DFS qubits. 



7. REDUCTION OF ALL REMAINING DECOHERENCE ON A SINGLE DFS 
QUBIT DURING LOGIC GATE EXECUTION 

The reduction of differential dephasing errors, as in the previous subsection, is particularly relevant for storage 
errors. However, this is only the first step. Additional sources of decoherence may take place during storage. 
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and in particular during the execution of logic gates. It is useful to provide a complete algebraic classification of 
the possible decoherence processes. This will allow us to see what can be done using decoupling pulses. To this 
end let us now write the system-bath Hamiltonian on two physical qubits in the general form 

HsB = -ffLeak + -f^Logi + -^DPS (39) 

where 

.ZI + IZ XY + YX XX -YY 
-Hdfs = bpanj , , 

-f^Leak = Span{X/, IX, YI, lY, XZ, ZX, YZ, ZY} 

= Spa.(X = ^^,y = ^M,^=^, ,40) 

where / is the identity operator, XZ = X1Z2 (etc.), and where Span means a linear combination of these 
operators tensorcd with bath operators. The 16 operators in Eq. (40) form a complete basis for all 2-qubit 
operators. This classification, first introduced in Ref. 16, has the following significance. The operators in i/oFS 
either vanish on the DFS, or arc proportional to identity on it. In cither case their effect is to generate an overall 
phase on the DFS, so they can be safely ignored from now on. The operators in -ffLcak arc the leakage errors: 
terms that cause transitions between states inside and oiitsidc of the DFS. A universal and efficient decoupling 
method for eliminating such errors, for arbitrary numbers of (c;iicodcd) qubits was given in Ref. 18. Finally, the 
operators in -ffLogi have the form of logic gates on the DFS. However, these are undesired logic operations, since 
they are coupled to the bath, and thus cause decoherence. 

In the previous subsection we showed how to eliminate the logical error Z, but we see now that this was only 
one error in a much larger set. To deal with the additional errors it is useful at this point to introduce a more 
compact notation for the pulse sequences. We denote by [r] a period of evolution under the free Hamiltonian, 
i.e., U{t) = exp{—iHsBT) = [t], and further denote 

P = i7i2(-|,0) = exp(-i|Xi2). (41) 

Thus Eq. (33) can be written as: 

exp[-z(i?J + BI){Z, + Z2)t] = [r, P, r, pt]. (42) 

As a first step in dealing with the additional errors, note that the symmetrization procedure [r, P, r, P^] can 
in fact achieve more than just the elimination of the differential dephasing Zi — Z2 term. Since X12 also 
anticommutcs with Y12 = \{YiX2 — X1Y2) G i?Logij if such a term appears in the system-bath interaction it too 
will be eliminated using the same procedure. 

So far we have used a ■^Xi2 pulse. Interestingly, the Hamiltonian X12 can also be used to eliminate all 
leakage errors. To see this, note that f7i2(±7r, 0) = exp(±i7rXi2) = Z\Z2. This operator anticommutes with 
all terms in i?Leak- Hence it too can be used in a parity-kick pulse sequence, that will eliminate all the leakage 
errors. 

At this point we are left with just a single error: X\2 ® B itself, in i?Logi- Clearly, we cannot use a pulse 
generated by X12 to eliminate this error. Instead, to deal with this error we need to introduce one more pulse 
pair that anticommutcs with X12, e.g., exp(±z-|yi2) = [/i2(±-|, ^). 

Let us now see how to combine all the decoherence elimination pulses into one efficient sequence. First we 
introduce the abbreviations 

H = I7i2(±7r, 0) = exp(±i7rXi2) = H^ = PP 

Q = t/l2(-|,|)=CXpH|Fi2) 

A = C7i2(±7r, |) = exp(±mFi2) = = QQ (43) 



9 



As argued above, the n pulse 11 eliminates Hi^eak- 



exp[-i(i?Logi + Hufs)2t] = [t, n, T, n]. (44) 

Now let us discuss adding the extra pulses needed to achieve full decoherence elimination. The 7r/2 pulse P 
eliminates Y and Z in ifLogi- Combining this with the sequence for leakage elimination we have the sequence of 
4 pulses: 

g-i(HDPS+^®Bjt)4r ^ [UiT)UU{T)Il]P^[U{T)UU{T)U]P 

= [r,n,T,P,T,n,r,pt], (45) 

(where we have used UP"^ = P, HP = P^). 

If we wish to entirely eliminate decoherence then we are left just with getting rid of the logical error due to 
X. To eliminate it we now combine with the F-direction, ■jt/2 pulse, Q: 



g-iHDPs8r ^ [U{t)ILU{t)PU{t)UU{t)P'']Q''[U{t)ILU{t)PU{t)UU{t)P'']Q 

= [r, n, T, P, T, n, r, Pt , Qt , n, r, P, r, n, T, Pt , Q] (46) 

which takes ten pulses. Unfortunately it is not possible to compress this further, since P^Q = {iX){—iY) = 
iZ and P^Q^ = —iZ, neither of which cannot be generated directly (in one step) from the available gate 
Uij {0, ) = cos 6*/ + i sin 0X^^.. . 

8. ELIMINATION OF OFF-RESONANT TRANSITIONS 

One important caveat in our discussion so far is that, because we need very strong and fast pulses, our gate 
operation may become imperfect. Specifically, off-resonant coupling may become important. This can cause 
unitary leakage errors from the DFS. These can in turn be reduced using the methods in Refs. 64, 65. Here we 
present an alternative and new method, that uses BB pulses generated in terms of the system Hamiltonian. 

Consider an A''-level system Hamiltonian 

JV 
i=l 

where the levels i = 1,2 denote our qubit. If all parameters are fixed, the evolution Uo{t) = exp{—iHot) is always 
on. If we turn on an interaction Hj, then the total Hamiltonian is H = Hq + Hj, which can be written as 

JV JV 

H = Y,Hu \i) {i\ + ^ {Hii \i) (il + h.c.) + Hl 

i=l i>j>2 

where the leakage term is 

JV N 

HL = \l)Y,Hu{i\ + \2)Y,H2il^\+ h.c. 

1=3 i=3 

Now note that 

JV JV 

U^mUoit) I*) + E {e-'^'^'-'^'^'Hij |z) {j\ + h.c). 

i=l i>j=l 

Using this we now show how to eliminate each of the leakage terms one by one. First, we eliminate |1) H13 {3\+h.c. 
by the BB sequence 

2H' = H + ul{-^)HUo{-^) 
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The new Hamiltonian H' does not contain |1) Hiz (3|+/i.c. Next, we eliminate |2) ifjs (3|+/i.c. by Uo{ ^^^^^ ), where 

^23 — t;xp(— i ^^'^J"^^^ )H23- We can clearly repeat this procedure so as to eliminate all leakage. The crucial 
point is that since we have used the system Hamiltonian to generate BB pulses, off-resonant transitions do not 
take place: Hq has no matrix elements between different levels. 

A concern is how we to obtain i7j(i) = exp{+iHot). This problem is shared by other methods dealing with 
the same problem. In principle it can be solved provided the level spacings are rationally related. Note in this 
context that usually we do not need to eliminate all Hu, since Hu ^ -ffii+i and typically decrease exponentially 
according to time-independent perturbation theory. We have assumed that there is no degeneracy. If there is, 
then all degenerate transitions will be eliminated simultaneously, so the procedure is simplified. 

9. COMBINING LOGIC GATES WITH DECOUPLING PULSES 

So far we have discussed computation using the encoded recoupling method (Section 4), and encoded decoupling 
(Sections 6,7). We now put the two together in order to obtain the full ERD picture. At least two methods 
are available for combining quantum computing operations with the sequences of decoupling pulses we have 
presented above. For a general analysis of this issue see Ref. 35. 

9.1. Fast + Strong Gates Method 

The decoupling pulse sequences given in Sec. 7 "stroboscopically" create collective dephasing conditions at the 
conclusion of each cycle. As noted above, this is equivalent to a periodic projection into the DFS. This property 
allows for "stroboscopic" quantum computation at the corresponding projection times.''''' Here the computation 
pulses need to be synchronized with the decoupling pulses, and inserted at the end of each cycle. Because 
of the conditions on the validity of the BB method, ^'^ the amount of time available for implementation of a 
logic gate is no more than the bath correlation time Tc = 2n/ujc. (An exception to this rule is the case of 1// 
noise. ^'') Assuming the dominant decoherence contributions not accounted for by the DFS encoding to come 
from differential dephasing (setting the Tc time-scale), and given that we already assumed that we can use pulses 
with interval At <^ Tc, it is consistent to assume that we can then also perform logic gates on the same time 
scale. 

9.2. Fast + Weak Gates Method 

There may be an advantage to using fast but weak pulses for the logic gates, while preserving the fast + strong 
property of the decoupling pulses. To see how to combine logic gates with decoupling in this case, let us denote 
by Hs = X^-X^. the controllable system Hamiltonian that generates the entangling gate Uij{9,4>i^4'j) [recall 
Eq. (15)]. Suppose first that we turn on this logic-gate generating Hamiltonian in a manner that is neither very 
strong nor very fast, so that the system-bath interaction is not negligible while Hs is on (this obviously puts 
less severe; demands on experimental implementation). Then the corresponding unitary operator describing the 
dynamics of system plus bath is: 

U{t) = exp[-zi(7Js + HsB + Hb)]. (47) 
Now, if we choose Hs so that it commutes with the decoupling pulses, then we can show that after decoupling 

U{t) ^ cxp[-i2t{Hs + Hb)], (48) 

provided t is sufficiently small. Tracing out the bath then leaves a purely unitary, decoherence-free evolu- 
tion on the system. To prove this, assume we have chosen t' and the decoupling Hamiltonian H'g so that (i) 
exp{-it'H's)HsBexp{it'H'g) = -Hsb, and (ii) [H'g,Hs] = . Then 

U{t)e-'''"sU{t)e-^''"'^ = c/(t)e-^*[^s+e-"'^sifsBe«'«fe+ifB] 

_ ^-it(Hs+HsB+HB)^-it(Hs-HsB+HB) 
^ ^-{2it{Hs+HB)+t'{[HsB,Hs] + [HsB,HB\)+0{t')} 

where we have used the Baker-Campbell-Hausdorff formula, exp(aA) exp(aS) = ex.p{a{A + B) + ^[A,B\ + 
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(49) 
(50) 
(51) 



Let us now show how to efficiently combine logic operations and decoupling pulses. For simpHcity consider 
only the case where we can neglect the X error, i.e., our decoupling sequence is the 4-pulse one given in Eq. (45). 
Suppose wc wish to implement a logical X operation, i.e., exp(— 1^X12). Recall that this involves turning on 

the Hamiltonian Hg = flxX^X^j, 17x^12 between two physical qubits. Because the decoupling pulses 
P = exp(— i|Xi2) and IT = exp(±i7rXi2) are generated in terms of the same Hamiltonian, they commute 
with Hg while eliminating Hsb (except for the terms in Hsb that have trivial action on the DFS). Thus the 
conditions under which Eq. (48) were shown to hold are satisfied. This allows us to insert the logic gates into 
the four free evolution periods involved in the pulse sequence of Eq. (45). Thus, the full pulse sequence that 
combines creation of collective dephasing conditions with execution of the logic gate is: 

g-it(nxXi2+HDPs) = Uit/A)TlU{t/A)PU(t/A)IlU{t/A)P^ , (52) 

with U{t) — exp[—it{Hg + Hsb + Hb)], and which, using the DFS encoding, is equivalent to the desired 
exp{— 10X12). This involves 8 control pulses, 4 of which are of the fast+strong type (those involving P and 11), 
and 4 of which must be fast, but need not be so strong that we can neglect Hsb- 

If we wish to implement logical Y operation, i.e., exp{—i9Yi2), then we cannot now use P and n, since they 
anticommute with Y12 and will eliminate it. Instead we should use decoupling pulses generated in terms of F12, 
which will also have the desired effect of eliminating Hi,eak, as well as X and Z logical errors, while commuting 
with the Y logic operations (and for this reason can of course not eliminate Y errors). These are just the Q and A 

pulses defined in Eq. (43). In gate terms this implies turning on the Hamiltonian H^ = f2YX^X0_|_^/2 ^yY 12 
between two physical qubits. Thus: 

g-it(n^ri2+HDFs) = u{t/4)AU{t/4)QUit/4)AU{t/4)Q\ (53) 

with U{t) = exp[—it{Hs + Hsb + Hb)], and which, using the DFS encoding, is equivalent to the desired 
exp(-i6'yi2). 

Finally, to generate single DFS-qubit rotations about an arbitrary axis we can combine Eqs. (52), (53) ac- 
cording to the Euler angles construction.^^ Given that Eqs. (52), (53) each take 8 pulses, the Euler angle method 
will generate an arbitrary DFS-qubit rotation in at most 24 pulses. 

Concerning gates that entangle two DFS-qubits, the situation may be more involved, since now the next- 
nearest neighbor pulses in Eq. (38), that create the collective dephasing conditions on four qubits, do not all 
commute with the U4 gate of Eq. (22). Therefore here we must resort to the strong + fast method of the previous 
subsection, i.e., we need to synchronize the U4 pulses with the end of the decoupling pulse sequence. We do not 
here analyze the situation with respect to the exchange Hamiltonian implementation of conditional logic gates. 

Taken together, the methods described in this section provide an explicit way to implement miivcrsal QC in 
a manner that offers protection against all sources of qubit decoherence, using a fast -|- strong (or fast -|- weak) 
versions of logic gates. 

10. DISCUSSION AND CONCLUSIONS 

We have proposed a method of encoded recoupling and decoupling (ERD) for performing decoherence-protected 
quantum computation. Our method combines an encoding into qubit-pair decoherence-free subspaces (each pair 
yielding one encoded qubit), with sequences of recoupling and decoupling pulses. The qubit encoding protects 
against collective dephasing processes, while the decoupling pulses symmetrize all other sources of decoherence 
into a collective dephasing interaction. The recoupling pulses are used to implement encoded quantum logic 
gates, either during or in between the decoupling pulses. 

The dynamical decoupling method requires an exponential number of pulses if the most general form of 
decoherence is to be suppressed, that can couple arbitrary numbers of qubits to the environment (total decoher- 
ence^^). This exponential scaling is avoided here by focusing on decoherence elimination inside blocks of finite 
size (e.g., at most four qubits) where arbitrary decoherence is allowed. However, we have implicitly assumed that 
there are no decoherence processes coupling different blocks. This is a reasonable assumption for most quantum 
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computer implementation, provided the different blocks can be kept sufficiently far apart until they need to be 
brought together in order to execute inter-block logic gates. When this happens, ERD can still be efficiently 
applied on the temporarily larger block. 

It may be questioned whether there is any advantage in using ERD compared to methods of active quan- 
tum error correcting codes (QECC). Both ERD and QECC are capable of dealing with arbitrary decoherence 
processes, and arc fully compatible with universal quantum computation. There are two main advantages to 
ERD: First, we need only two physical qubits per logical qubit, compared to a redundancy of five to one, in 
order to handle all singlc-qubit errors in QECC.^ So far the most advanced experiments outside NMR, i.e., those 
involving trapped ions, have used up to four qubits,®" so that this encoding economy is a distinct advantage for 
near-term experiments. Second, our method is directly compatible with Hamiltonians describing a variety of 
quantum computer proposals. On the other hand it is not clear how to directly use QECC given Hamiltonians 
describing specific systems. These are general features of ERD: economy of encoding redundancy and use of only 
the most easily controllable interactions. The disadvantage of ERD compared to QECC is that there does not 
exist, at this point, a result analogous to the threshold theorem of fault tolerant quantum error correction. This 
means that wo cannot yet guarantee full scalability of ERD as a stand-alone method, because wc do not yet 
know how to compensate for imperfect pulses. However, in principle it is always possible to concatenate ERD 
with QECC, as done, e.g., for DFS with QECC in Refs. 66, and then the standard fault tolerance results apply. 
We expect that the theory of composite pulses®'' will also play a key role in this further development of ERD. 

We hope that the methods proposed here will inspire experimentalists to implement encoded recoupling and 
decoupling in the lab, thus demonstrating the possibility of fully decoherence-protected quantum computation. 
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